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NOTATION 

Function  representing  the  ship  surface 
Function  representing  the  free  surface 
Acceleration  of  gravity 
Defined  in  C  73  ^  C9l>  and  [ill 

Green's  function 
Draft  of  the  ship 

=  g/v3 

Length  of  the  ship 

Intersection  of  the  plane  C  =  0  and  r\  =  f(?,C) 
Normal  vector  on  the  surface  into  the  flow 
Pressure 

Defined  in  [8],  [ 10] ,  and  [ 12] 

Wave  resistance 
Free  srrface 

Ship  surface 

Uniform  velocity  at  infinity 

Right  handed  rectangular  cartesian  coordinates 
as  shown  in  Figure  1 

5  (?  -  )  cos  9  +  (t)  -  y1  )  sin  6 

Fictitious  friction  force 

Small  parameter  representing  the  beam  length  rati 
of  a  ship 

Transformed  coordinate  system 
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represents  coefficients  of  In  the  expansion 
of  the  attached  quantity 


represents  the  corresponding  quantity  in  the 
transformed  space 
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SUMMARY 

The  free  surface  condition  of  the  surface  wave  is  con¬ 
sidered  on  the  free  surface  itself  instead  of  the  mean  free 
surface,  by  the  use  of  a  coordinate  transformation,  together 
with  the  scheme  of  a  systematic  expansion  in  a  small  parameter. 
Thus  a  higher  order  ship  wave  theory  is  developed.  The  most 
common  practical  case  of  a  slender  ship  with  an  almost  flat 
bottom  is  especially  treated  in  detail.  The  lowest  order  re¬ 
sult  is  the  same  as  that  given  by  Michell's  theory.  The  next 
higher  order  potential  and  wave  resistance  are  derived  ex¬ 
plicitly. 
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INTRODUCTION 

Since  Kelvin  (1887)  found  the  solution  for  linear  water 
surface  waves  and  Michell  (1898)  formulated  the  wave  resistance 
due  to  a  thin  surface  ship  advancing  in  an  inviscld  fluid,  many 
experimental  and  theoretical  works  on  this  subject  have  been 
performed  by  various  ship  hydrodynamists .  Th  waves  and  the 
wave  resistance  of  many  ships  have  been  calculated  and  compared 
with  experiments.  The  theory  of  minimum  wave  resistance  has  been 
developed  and  waves  from  the  bow,  the  stern  and  the  shoulder  of  a 
ship  have  been  analyzed. 

These  theories  are  based  on  solutions  of  the  Laplace  equa¬ 
tion  with  linearized  boundary  conditions  on  both  the  free  sur¬ 
face  and  the  ship  surface.  The  theories  are  mathematically  ele¬ 
gant  and  certainly  very  significant  in  the  development  of  the 
ship  wave  theory.  However,  the  general  agreement  with  experiment 
has  been  relatively  very  poor  (especially  for  practical  ships), 
as  compared  with  other  linear  theories  such  as  have  been  devel¬ 
oped  in  aerodynamics  and  for  cavity  flows. 

To  improve  this  situation,  the  development  of  higher  order 
wave  theory  has  been  suggested  in  terms  of  a  systematic  expan¬ 
sion  in  small  parameters  (Stoker,  1957).  Sisov  (1961)  actually 
formulated  the  second  order  wave  resistance  for  a  thin  ship  where 
the  beam-length  ratir  '  >  considered  to  be  a  fundamental  small  pa¬ 
rameter  as  in  Michei*  -o  theory.  To  improve  the  comparison  be¬ 
tween  wave  theory  and  experiment,  the  streamline  tracing  tech¬ 
nique  was  recommended  (inul,  1957),  and  also  a  slender  ship 
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theory  has  been  developed  (Vossers,  1962;  Tuck,  1963;  Maruo, 

1962;  Joosen,  1964).  Wehausen  (1964)  has  also  considered  an  ex¬ 
act  formulation  of  the  ship  wave  problem  and  showed  that  for  a 
thin  ship  with  small  draft,  the  most  important  higher  order  ef¬ 
fect  can  be  represented  by  a  line  integral  along  the  intersec¬ 
tion  of  the  ship  surface  and  the  water  surface. 

Although  no  proof  has  been  made  for  the  convergence  of  the 
expansion  of  the  potential  in  power  series  of  a  small  parameter 
as  in  Sisov  or  Wehausen' s  work,  this  expansion  scheme  is  common 
in  applied  mathematics.  (Lighthill,  1934,  Wehausen,  i960. 

Van  Dyke,  1965).  The  usual  expansion  problems  are  made  more  dif¬ 
ficult  in  the  case  of  ship  waves  because  of  the  unknown  position 
of  the  free  surface  boundary  where  it  is  necessary  to  apply  the 
free  surface  boundary  condition.  In  addition,  the  free  surface 
is  internally  bounded  by  the  ship  surface  where  the  usual  rigid 
boundary  condition  must  be  applied. 

To  resolve  the  difficulty  connected  with  the  unknown  free  sur¬ 
face  position,  here  we  first  transform  the  coordinates  such  that 
C  =  0  always  represents  the  free  surface.  Then,  we  apply  the 
scheme  of  an  expansion  in  small  parameters  to  all  of  the  physi¬ 
cal  quantities;  in  this  development,  first  the  beam-length  ratio 
was  assumed  small  and  then,  in  addition,  the  draft-length  rati  , 
was  assumed  small  (i.e.,  the  ship  Is  thin  and  then  slender). 

After  substituting  the  relevant  series  in  our  governing  equation 
and  the  boundary  conditions,  we  equate  all  the  same  order  terms 
in  each  equation.  Then  the  governing  equation  which  is  the 
Laplace  equation  In  the  physical  space  becomes,  in  the  transformed 
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space,  Laplace  equations  in  the  lower  orders,  and  Poisson’s 
equations  in  the  higher  orders.  In  each  order,  the  problem  is 
a  well  defined  linear  boundary  value  problem  whose  solution  is 
available  by  the  use  of  Green’s  theorem,  with  one  well  known 
Green’s  function  which  is  equivalent  to  a  potential  due  to  a 
source  under  the  linearized  free  surface.  The  higher  order  so¬ 
lutions  depend  upon  the  lower  order  solutions.  Therefore,  we 
first  solve  the  lowest  order  problem,  and  use  this  solution  for 
the  next  higher  order  problem,  and  so  on.  Thus,  the  potential 
and  the  wave  resistance  for  both  thin  ships  and  slender  ships 
are  formulated  explicitly  here  in  both  the  first  and  the  next 
higher  order. 

The  lowest  order  solution  is  the  same  as  Michell’s  solu¬ 
tion  whether  the  ship  Is  thin  or  slender,  as  Wehausen  showed. 

The  next  higher  order  solution  is  relatively  simple  especially 
for  the  case  of  fast  slender  ships.  This  is  due  to  three  ef¬ 
fects:  the  effect  from  the  line  integral  (earlier  mentioned) 
which  Is  due  to  tne  free  surface  bounded  internally  by  the  ship; 
the  effect  from  the  change  of  submergence  due  to  the  wave;  and 
the  bottom  effect.  The  second  order  solutions  for  the  cases  of 
thin  ships,  and  of  slow  and  medium  speed  slender  ships  (ora inary 
merchant  ships)  are  a  little  more  complicated  since  the  non¬ 
linear  effect  cf  free  surface  needs  to  be  included.  Singulari¬ 
ties  at  the  bow  and  stem,  especially  at  the  intersection  with 
the  free  surf  ice  are  carefully  analyzed.  For  simplicity,  no 
trim  nor  sir  age  is  considered,  although  these  effects  can  be 
added  witho  t  too  much  difficulty. 


w  auk*. 
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SISOV’S  HIGHER  ORDER  WAVE  THEORY 

We  consider  the  rleht-handed  rectangular  coordinate 
0-x,y,z  or  0-xx  ,yx  ,zx  and  the  surface  ship  y  =  f(x,z)  as  in  Fig¬ 
ure  1,  in  a  homogeneous,  inviscid,  Irrotational,  and  infinitely 
deep  fluid  with  a  free  surface.  The  flow  at  infintiy  is  con¬ 
sidered  to  be  uniform  with  the  velocity  U,  and  there  exists  a 
perturbation  velocity  potential  cp,  which  satisfies 

Yacp  =  0  Cl] 

and  the  proper  boundary  conditions  on  the  free  surface  and  the 
ship  surface,  with  ?cp  =  0  and  3p  =  C  at  Infinity.  Then  by  Green's 
formula 


<!P 


_1_ 

4tt 


*(xi  ,yx 


t  <-i 


)  Gn(xx  ,yi 


;x,y,z) 


-  <P_(xi  fYi  ,.Zi  )  G(xi  ,y: 

f  a 


,y,z)  dS 


where  the  subscript  n  indicates  a  partial  derivative  in  tr.e  di¬ 
rection  of  the  inner  normal  Into  the  fluid  at  the  whole  bound¬ 
ary  S(xx,yx,2x)  including  the  free  surface  Sv  and  the  ship  sur¬ 


face  3  in  Figure  1.  0  is  the  Green's  function,  which  is  bar 

s 

mor.ic  everywhere  in  the  fluid  except  at  *  point  (x,v.u)*  (x*  ,y- 
where  it  has  a  singularity  lixe 


■*  ) 
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and 


l/  {(X-X;)8  +  (y-JTi  )’  +  (Z-Zi)8} 

> 

VG  =  0  ,  and  G  =  0  at  infinity 

in  z  <  0 


The  boundary  conditions  for  cp  on  the  free  surface  z  =  F(x,y) 
are:  the  kinematic  condition 


(V-cp  )  F  -  cp  F  +cp  =  0 
x'  x  y  y  z 


[4a] 


and  the  dynamic  condition 


F(x’y)  -  i  ’x  +-k  (<P*  V*"**  > 


=  0 


[4b] 


The  boundary  condition  on  the  ship  surface  y  =  f(x,z)  i. 
the  kinematic  condition 


or 


(V-cp  )f  +  q>  -  cp  f  =  C 
x '  x  y  '  z  z 


cp  -  n  •  V  =  0 
n 


[5] 


For  convenience  we  do  not  consider  here  any  trim  or  sinkage  of 
the  ship.  We  consider  a  thin  ship  symmetric  with  respect  to  the 
v  =  0  plane,  with  small  slopes  to  the  y  -  0  plane,  and  a  small 
value  of  the  half  beam-length  ratio,  e;  and  assume  the  series, 
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cp(x,y,z,e)  =  ecPi  (x,y,z)  +  eacpa  (x,y,z)  +  — 

P(x,y,c)  =  epx  (x,y)  +  e3Fa  (x,y)  +  — 

f(x,z)  =  efi  (x,z)  [6] 

By  substitution  of  [6]  in  [4a],  [4b],  and  [5],  and  using 
Taylor's  expansion,  we  obtain  (see  e.g.  Sisov,  1961)  comparing 
the  terms  of  e^, 

+tpJxx  =  pj(x’y)  0nSPF  C7] 

(kQ  =  g/v3 ) 

for  the  free  surface  condition 
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Pa  (x,y)  =  —  ^plxa  +  cpiya+  ~  cpizaj 


-  I  ** 


|kQCPiZ  +  CPl  XX 


■) 


Cn] 


q a  (x,z) 


[12] 


In  general. 


Pj(x,y)  and  q^(x,z)  for  J  * 


2  are  represented  In 


terms  of  known  functions  Including  the  obtained  lower  order  solu¬ 
tions  {^cpi  ;  i  <  j}  . 


Sisov  (1961)  obtained  two  Green's  functions  for  the  free 
surface  singularities  on  z  =  0  and  for  the  ship  surface  singu¬ 
larities  on  y  =  0  respectively,  both  satisfying  the  first  order 
free  surface  condition  [7]  and  [9]  except  at  the  singularities. 
Thus,  if  we  consider  that  the  boundary  conditions  everywhere 
on  the  z  =  0  plane  and  the  y  =  0,  z  ^  0  haif  plane  are  given  as 
in  [7]  and  [8],  respectively,  the  combination  of  the  two  Green's 
functions  will  lead  us  to  the  solution  of  Sisov  (1961). 


The  problem,  however,  is  that  the  free  surface  is  only  out¬ 
side  of  the  ship  boundary  and  the  actual  submergence  of  the  ship 
depends  upon  the  wave  height.  Thus,  the  boundary  conditions  on 
Spp  and  S  s  are  not  necessarily  those  given  everywhere  on  the 

z  =  0  plane  and  on  the  half  y  =  0  (z  *  0)  plane.  In  addition. 

In  practice,  our  ships  are  not  thin  and  the  slope  to  the  x,z  plan 
is  not  everywhere  small. 
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To  overcome  these  difficulties  we  will  try  another  formu¬ 
lation  of  the  same  problem  with  a  different  approach. 

COORDINATES  DEFORMATION  AND  FORMULATION  OF  PROBLEM 

We  consider  a  transformation  of  the  (x,y,z)  coordinates 
to  the  (§,ti,C)  coordinates  by 


x  =  § 
y  =  r\ 

z  =  C  +  F(x,y) 


[  13^ 


Thus,  C  =  0  represents  the  free  surface  z  =  F(x,y).  Of  course, 
at  the  beginning,  we  do  not  know  the  wave  height  F(x,y).  Our 
governing  equations  and  the  boundary  conditions  [  l] ,  [4]  and  [5] 
can  be  represented  in  terms  of  the  ?,T|,C  coordinates  through 


substitution  of 


cp  (x,y,  z )  =  cp{§,T),C  +  F(x,y )}  =  cp#(?.R;C) 

’x  =  q,g  -  -  *t*X 

*xx  “  *55  -  ^c'x  -  *CFxx  + 


cp  =  cp*  -  cp*F 

y  u  C  y 

cp  =  cp*  2cp*,  F 
yy  r\r\  y 


-  cp*F  +  cp*  F  8 

C  yy  CC  y 


[14] 


cp  =  cp* 

2  C 


*zz  =  VCC 
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in  Cl],  [4],  and  [5].  Now  if  we  use  the  series  [6]  in  the 
transformed  space 


Acpj  = 

ko*jC  +  I’jSS  =  PJ  (5<tl) 

on 

S*F 

> 

[151 

on 

S* 

S 

J 

where  S^*  and  S*  indicate  the  corresponding  surfaces  of  and 

S  in  the  transformed  soace . 
s 


p*  =  0  ) 

q*  =  -Vf* 

1  1  ?  J 


[16] 


f*(5>C)  s  f*Cx,z  -  F(x,y)}  *  f(x,z) 

F(?,r|)  =  F(x,y)  -  |<px  -  ^  {cpxa  +  cpya  +  V*  )  l 
Fi  (?.l)  =  |  P*  on  C  -  0 


[  17  ] 
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i?rix  T  “"lnTiy  T  '•'lC  (Flxx  +  Flyy) 


da  (S,ri,C )  =  ap%,Pn„  +  a»*  P,„  +  cp*  b .  +  P. 


P?(§,n)  =  A  +  cp«  +  i cp*’  ]  +  (ap*?cF 


+  CD*  F 
lx  1C  Ixx 


\  L l3al 


#  / 

Qa  (§,C  )  =  n P*-f*.  +  cp*  f*  ^  +  Vf*  F  +  cp*  F 
v  ,w  ^l?  l?  ^ic  1C )  1C  lx  vic  ly 


Fa  (x,y )  = 


% 

g  ^2x 


-if- 


3  +  cp  3 
lx  *ly 


cp 


Iz 


V  * 
—  cp* 


g 


F 

1C  lx 


+  cp 


*3 

m 


<) 


[  l8b] 


V  pr 


for  j  ^  2  are  represented  In  terms  of  known  func¬ 


tions  Including  the  obtained  lower  order  solutions,  {cp*  ;  1  <  j}. 
We  drop  the  superscript  *  except  In  the  case  of  confusion. 


For  a  general  ship  which  has  a  small  draft-length  ratio  as 
well  as  a  small  beam-length  ratio,  the  boundary  condition  cor¬ 
responding  to  [8]  and  the  order  analysis  needs  to  be  dealt  with 
carefully.  The  exact  boundary  condition  on  the  ship  surface  In 
x,y,z  coordinates  is 


Vf 

_ x _ 

V  i+f 3  +f 3 

V  x  z 


on 


S 

s 


[19] 
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By  transformation  to  the  5,rbC  coordinates 


♦v  = 


V 


1+ 


l+f,3+f  8 

—2 - a -  ^Cp.F  -  F  fcp.f»  +  cp-f 

a  a  ]C  y  x  C  5  '  C 


f.-f-F  ]s 
15  C  x 


+f, 


Vf.+Vf-F  >  on  S* 

5  C  xf  s 


5 


[20] 


SOLUTION  FOR  EACH  ORDER  POTENTIAL 

The  solution  for  each  j  can  be  expressed  using  Green’s 
formula  in  the  transformed  space. 


*  Tj^p  J  J  Uj(xi  ,yi  ,zi  )  Gv(xi  ,yi  ,Zi  ;5,n,C 

s*  L 


-  v (xx  ,yi  ,zi  )  G(xi  ,yi  ,zx  ;§,h,C  )  J 


dS 


-  h  JJJ 


G  AcPj  dt 


D* 


where  the  Green's  function  G  satlsfie; 


[21] 


G^  +  *0  on  C  *  0 


[22] 


0  on  r\  m  o 
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in  addition  to  the  usual  properties  previously  mentioned  in  [3] . 
Such  a  Green's  function  G(xi  ,yx ,  zx  j?  ,ruC  )  is  w*»ll  known  to 
represent  the  potential  due  to  a  point  source  located  at 
(xi,yi,zi),  (see  e.g.,  Lunde,  1952). 


G(xi  ,yi  ,zi  ) 


J_  _  _1 
n  ”  ra 


7T  * 


-^Reff  2«£»l 

r  J  J  k-k  S( 


-k|  C+Zi  |  +  ikuu 


seca0  -jji  seed 


dkd0 


-7T  o 


_1_ 

ri 


i 

ra 


t/2  -koseca0{  |  C+zi  | } 
+  4k  f”  e 

-tt/2+6 


sec8 0 sin(kotuseca  0  )d0 


-ir/2+6  0 


mu,tkoseca 0sin(m|C +zx  |  -mcos(m|C+zxl )} 

k  8  sec4  0  +  ma 
o 


[23] 


where 

(?-Xx)a  (*n-yr  )a  +  (C-Zi)a 
r*a  =  (5-xx  )a  +  (r|-y,  )a  +  (C  Zl  )a 


6  «  arctan 

«  (5-Xx  )  COS0  +  (ti-y^  )  sln0 

M  a  fictitious  friction  force  which  is  put  to  zero  after 
integration . 
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G  (xi  ,yx  ,0;?  ,ti,C  )  represents  the  potential  due  to  a  Kelvin’s 

Xl 

source  or  a  pressure  point  on  the  free  surface.  Using  the  con 
dltions  [22 ]  for  Green's  function  and  the  boundary  conditions 
for  we  obtain 


S* 

F 


Vv'V3  dS  ■  ■  Va  •  ^  dXldyi 


fair  VW  +  rp?  dXldyi 


s* 

F 


7—  cp.G 

ko  J  Xl 


’j*°  dyi  +  * 


p*G  d  >'i  ciy; 


[  23*1 


S*F 


where  l  represents  the  intersection  of  S*  and  .  Therefcre  we 
can  write. 


cp  =  E  e°cp,  ^ 


JL  £  g--1 

UTT 


s* 

s 


*.G  -  q*.  G  dS 
J  v  jv 


+  TT  €  0 


J  Xi  J  X;  dXx 


dfi  .  »  GciX;  iy, 

— *-  dXi  +  p*  - 


J 


S* 

F 


MjdT 


r  ..  -5 

L  *  J 


D* 
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Now  we  may  consider  a  Taylor's  expansion  of  the  Green's 
function  at  t|  ■  0, 


-a  f 

G(xi  ,efx  ,zi  )  «  G(xx  ,0,Zx  j5 G, 


Yiyi 


+  — 


G^(xi  ,efx  ,zi  ;5,h,C  ) 


-G  f  +  0  -  0  f 

-j&.  a _ & _ ^  z.i. 


a  +  f  a 

Zx 


>  [25] 


*  efiGv  „  (xx  ,0,Zx  ;?,T1,C)  +  0(ea),etc. 
yiyi 


j 


Since  we  do  not  know  a  priori  we  consider  the  sequence 
of  ship  surfaces  which  converges  to  Sj,  or 


s*08  •  9*  — 1  • 


[26] 


where  S£s  represents  the  ship  surface  without  any  wave.  Thus 
for  each  tp  we  can  have  a  sequence  of  solutions  which  may  even- 

v 

tually  have  the  domain  of  integration  S*.  The  problem  for  the 

3 

first  order  potential  on  S*  is  exactly  the  same  as  the  first 

os 


order  problem  in  the  case  of  Sisov's  theory. 
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SLENDER  SHIPS  WITH  AlflOST  FLAT  BOTTOM 


When  we  consider  that  the  draft  of  the  ship  is  of  the  same 
order  as  the  beam  in  [24],  the  draftwise  integration  of  the  in¬ 
tegrand  O(e^)  becomes  0(e*^+1).  Thus,  in  [24],  the  lowest  order 
of  <p  is  0(ca  )  on  and  outside  of  the  ship.  If  the  ship  has  al¬ 
most  a  flat  bottom,  q>  is  0(c)  only  near  the  ship  sides  (from 

C  20]  ) . 


With  all  these  considerations  including  Equations  [24]  - 
[26],  if  we  collect  the  coefficients  of  the  same  order  from 
[24]  we  can  have  for  each  j.  Thus  the  lowest  order  solution 
is 


L  -Hi 

~tf/ /f'*> 


G(xx,0,ezx  ;x,y,z)  dzxdxi 


0  0 


[2?] 


where  H  *  eHx  and  L  indicates  t.*e  draft  and  the  length  of  the 
sMp  respectively. 

The  next  higher  order  solution  would  be 


q>  •ca<pa-K3<pa  m  1- J  J  f?Xl  (xx  ,zx  )c(xx  ,0,zx  ;x,y  ,z-Fx  (x,y  )}dzx  dxx 

Sfpr\s 

£»  r  r  L 

*  2?  I  J  (xi  (xx  ,0,-H)C  (xx  *0,-H;x,y,2)  dxx 

*  ST  ,f:  '0)°Xl  .O.y.y.z).  ^»4*,  .fi  .C)gJ|£-  ^x,J 


Ue] 
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The  first  term  of  the  right  hand  side  in  [28]  is  different  from 
[27]  only  In  the  domain  of  integration.  Thus,  the  influence  of 
the  change  of  the  submerged  ship  surface  due  to  the  wave  ap¬ 
pears  in  the  change  of  flx(x,z)  on  S*„  into  fi$(S,C)  on  Sfs 

which  is  bounded  by  a  straight  plane  C  =  0  at  the  exact  free 
surface  and  the  wavy  ship  bottom  in  the  §,h»C  space.  If  the 
ship  is  wall  sided,  the  only  influence  is  from  the  vicinity  cf 
the  bottom  but  not  from  the  free  surface  in  Sfs. 


The  second  term  in  the  right  hand  side  of  [28]  is  the  flat 
bottom  effect.  Here  the  Taylor  expansion  of  <f>  which  is  symmetric 
with  respect  to  h  *  0  Is  used.  It  represents  a  distribution  of 
vertical  doublets  with  strength  proportional  to  fjCPa  . 

The  third  integral  in  the  right  hand  side  of  [28]  is  the 
line  integral  along  the  line  of  intersection  between  the  snip 
surface  and  the  free  surface.  This  indicates  the  line  distrlt 
tlon  of  singularities  on  the  free  surface,  and  will  be  investi¬ 
gated  further  in  the  following  section.  This  is  particularly 
important  because  even  for  a  thin  ship  with  deep  draft,  cp  will 
be  of  the  same  ord«v  as  the  case  of  the  small  draft  when  the 
ship  speed  is  low. 


The  wave  height 


and  from 


1 J 


sufficiently  apart 
It  may  be  written 


fa  (  »y )  *  ~5ax(x»y>^) 
£ 


from  the  slender  ship.  On  the  ship  surface 
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F«(x,fi)  =  ifvtpaxfx.fi.O)  -  'y 


[29b] 


except  near  the  bow  and  the  stern. 


THE  LINE  INTEGRAL  ON  THE  FREE'  SURFACE 


The  third  integral  in  the  right  hand  side  of  Equation  [28] 


I  = 


CPs  (xi  ,fi  ,0)GXi  (xi  ,fi  ,0,?  ,P,C  )-  CP2X1  G  |  dxx  [30a] 


can  be  interpreted  as  a  potential  due  to  the  distribution  of 
doublets  and  sources  on  a  line  p  =  fj.(§,0)  whose  strengths  are 
proportional  to  cp2  fx  l/kQ  and  cp8  Xi  flx  x  respectively .  It  is 
well  known  that  a  doublet  distribution  on  the  free  surface  is 
the  same  as  a  pressure  distribution  on  the  free  surface  (Wehausen, 
1959;  Ursell,  i960).  Therefore  a  distribution  of  sources  on 
tne  free  surface  can  also  be  Interpreted  in  terms  of  a  corres¬ 
ponding  distribution  of  pressure  on  tne  free  surface.  For  ex¬ 
ample,  we  consider  a  smooth  distribution  of  pressure  p(S,p)  on 
a  certain  domain  D  on  the  free  surface  C  =  0  of  an  otherwise 
uniform  flow,  and  p  =  0  on  the  boundary  of  D.  Then 


V  = 


V 


•  tTTf  g 


pG  dxx  dyi  =■ 
Xi 


P  G  dxx  dyi 

Xi 
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Namely,  the  potential  anywhere  in  the  fluid  can  be  expressed 
either  by  a  doublet  distribution  on  £  =  0  with  strength  pro¬ 
portional  to  p(?,r|)  or  a  source  distribution  on  C  =  0  with 
strength  proportional  to  p^(§,T|).  Singularities  on  the  free 
surface  need  special  attention  because  it  is  known  that  when 
a  pressure  point  or  a  point  doublet  is  located  on  a  free  sur¬ 
face,  the  flow  behavior  is  singular  not  only  at  the  point  it¬ 
self  but  also  on  its  path  line  and  the  corresponding  wave 
resistance  blows  up  (Lamb,  1932;  Ursell,  i960).  The  same  kind 
of  singular  behavior  takes  place  for  a  point  source  on  the  free 
surface  and  the  corresponding  wave  resistance  becomes  also  in¬ 
finity  (Yim,  1966). 

By  an  integration  of  [30a]  by  parts 


I 


|L 


~  ^2  fix  Gr(xi  ,  fi  ,  0,§  ,  T],C  ) 
0 


Xi  =0 


[30b] 


neglecting  higher  order  terms.  The  first  term  of  the  right 
hand  side  in  [30b]  indicates  the  potential  due  to  a  point  source 
located  at  (xi,fi,0)  on  the  free  surface.  This  would  induce  the 
singular  behavior  mentioned  before  unless  T\xi  =  0  at  Xi  =0 
and  Xi  =  L.  Therefore,  in  this  report  we  will  consider  the  ship 
surface  y  =  f(x,z)  which  is  sufficiently  smooth  in  -®  <  x  <  ®, 
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near  the  free  surface.  In  an  engineering  approximation  of  the 
line  Integral  for  the  lower  orders,  the  free  surface  water¬ 
line  of  a  ship  of  finite  entrance  angle  may  be  approximated  by 
a  cusped  smooth  waterline  sufficiently  close  to  the  original 
one,  since  there  can  never  exist  such  singular  phenomena  in  the 
actual  case  of  real  fluid.  At  present,  no  meaningful  mathe¬ 
matical  limit  of  the  wave  resistance  or  the  wave  height  due  to 
either  a  point  source  or  a  point  doublet  on  the  free  surface  is 
available.  This  is  discussed  more  fully  later  in  the  section 
entitled  Singularities  at  the  Bow  and  Stern. 

WAVE  RESISTANCE 

The  wave  resistance  R  can  be  calculated  by  the  integration 
of  pressure  p  from  the  Bernoulli  equation, 

P  =  P  jv<Px  -  zg  -  |  (cpx3+  cpy2+  cpz2  )J  [31] 

over  the  wetted  ship  surface. 


o  -H 
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If  we  put 


R  =  e2  Ra  +  £3  R3  +  £4  R4  +£6R5  + - 


the  lowest  order  wave  resistance  for  a  thin  ship  will  be 


since 


Rs  = 


L  ^-o 

2pVcplxfixdx  dz 


This  leads  to  the  Michell's  wave  resistance  formula.  The 
order  wave  resistance  due  to  a  slender  ship  is 


R4 


2p(  Vcpa  x- 


fj.  xdx  dz 


in  which  only  the  part  due  to  regular  waves  will  remain, 
this  is  also  the  same  as  the  Michell's  wave  resistance. 


L  33  J 


C34] 


lowest 


[33] 

Thus 


1 
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The  next  higher  order  wave  resistance  of  a  slender  ship 
with  the  almost  flat  bottom  is 

L  o  L  Fa(x,z) 

Be  =  2  pVcp3Xfixdx  dz  +  2  /  I  pVcp2Xfixdx  dz 

o  -Hi  oo 


.L^Fa  (x,y) 


-2  /  p  zg  +  |  Vs  x2 


fi  xdx  dz 


[36] 


o  o 


where  cpa  ,  cp3  ,  and  Fa  are  given  in  [27],  [28],  and  [29b] 
respectively.  If  we  consider  that  fix  is  almost  independent  of 
z  near  the  free  surface,  the  last  integral  will  become 


j"  P  (Fsg  +  Vs  fiX8  ]  FafiXdx 


From  [28],  [35]  and  [36],  we  can  write 
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P  ■=  e4R4  +  ee  Re 


Vf?§d5dC 


S* 

Iphs 


f? Xi  G-  J  y.\  ,  0,  Z\  ;?jO,  «l  dXj.  dzi 


SlpP£ 


opys 


fixdxdz  fi  (xi  ,  -ri)cp3  (x-i  ,  0,  -H)G  (xj.  ,  0,  -H;x,  0,  z  )cix 

dV  Zix 


+  “  G^  Jcpa  (x-x  »  0,  0  )  fj.  xx  Xx  +  2%  (*i  ,  0,  0 )  fi  xi 


(xi,0,0)fi  XlJdxi 


„  L 

J  P  (Fsg  +  V*flx3  j  F2flxdx 


tv?] 


Since  f.x  =  0  or;  me  bott  om,  the  integration  over  the  com  or 
is  omitted .  However,  In  actual  calculation  it  may  be  letter 
corsldei  the  forma  1  consideration  of  the  bettors  integral  to 
e'Unnat  e  the  calculation  of  local  aisturoat.ee,  because  If  w-- 
or-  3  icier  the  wave  resistance  due  to  a  singularity  dlstr ifcut  lor. 

•  mo  .igr  I  aga  Li  ey  c  t  heorern  t  he  c  out  id  burl  or.  i  i  or  *  no  .  oc  a  . 
a  1  ::r  nrbat.ce  is  zero  (see  e.g.,  Lunde,  19V7)-  Ev  Up  same  reas  m 


»...v 
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we  consider  the  effect  of  cp2  x  on  the  line  source  distribution 
on  the  free  surface 


_2p  f  ~a~is-^--0-0)  j4a(x,0,0)  fixx  +  %x(x>0,0)fiX'jdx 


H 


■  -2p  I  7 


Pag  +  \  ^fix2  {  f2-  fixx  +  V 


Fa  +  pr-  fi  x3  [  x^dx 
o  1  j 


[38] 

If  we  add  [38]  in  [ 37 3  the  effect  of  the  local  disturbance 
due  to  the  surface  line  integral  is  eliminated  within  oCe6  ), 
and  subtracting  [38]  from  the  last  integral  of  [ 37l  we  obtain, 


r = *{r^- 


+  01 8  )  sec3 9  d9 


rvr/2 

-  2  e4  (C  C.  +SSJ  see6 9  de 

J  v  o  b  o  b' 

o 

,v/2 

-  2  I  €4  (C  C_  +SS.)  sec3  6  d9  \ 

J  O  F  c  !•  j 

o  1 

+  e4  n  1  f  cu  £  +  V  f  91  lfa  f  al  c  « 

+  e  p  j  '  “  1 1  xx  x:  *  a  7?  +  2  *x  !+  l,ag4-v  Ux 


\ 


— —  f> 


x6  ^  d  x 


[59] 
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where 


o  J  S 


■u 


dx  dzVk  fix  e 
o  A 


opys 


k  zsec80  r  ^ 
o  cos 


.sin. 


(k  xsecQ  ) 
v  o  ' 


Ci 

Si 


-Jf  d5dCVkof1?e 


k  c  sec®  6 

o 


St  w 

lpTlS 


COS 

sin 


(k  ?  sec0 ) 
o 


rL 

kQafi  (x,-H)cpa  (x,0,-H)  e 


-k  Hseca0 


o 


cos 


sin 


(k  xsec0  )d 
o 


V  L 


1 


F  J  o 


=  J  |<Pa  (x,0,0)fiXX+  2cpaxfiXJ  |  j(kQxsec0)dx 


Thus  the  higher  order  wave  resistance  of  slender  ships 
0(c5  )  is  from  the  change  of  the  ship  submergence  due  to  the 
wave,  the  bottom  effect,  and  the  line  integral  due  to  the  free 
surface  bounded  internally  by  the  ship  surface.  These  higher 
order  effects  are  functions  of  the  lowest  order  wave  height  i  nd 
the  lowest  order  potential,  which  are  zero  at  infinity. 
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ORDER  ANALYSIS  AND  DEPENDENCE  ON  k 

o 

The  selection  of  terms  In  each  higher  order  for  thin  ships 
and  slender  ships  in  the  above  analysis  has  been  strictly  for¬ 
mal  under  the  assumptions  that  series  expansions  [6]  converge 
for  any  small  e,  and  all  the  variables  behave  properly.  How¬ 
ever  this  analysis  would  be  incomplete  if  we  did  not  consider 
the  actual  situation  which  is  highly  dependent  on  the  Green's 
function  [23]  and  the  parameter  kQ.  In  practice,  for  ordinary 
merchant  ships  (low  and  medium  speed  ships)  and  fast  ships,  the 
magnitude  of  the  following  parameters,  viz,  Froude  number 
Fr  =  Wy  ,  ^qL0  s  1  beam-length  ratio,  and  draft -length 
ratio  are  as  shown  in  Table  1. 


TABLE  1 

Magnitude  of  Ship  Parameters 


Ordinary 
Merchant  Ships 

- 1 

Fast  Ships  ! 

F 

r 

0.33  —  0.18 

| 

0.33  —  1.5  j 

l/kQL  *  lAi 

0.11  —  0.033 

0.11  —  2.2  ! 

t 

k0L  *  kl 

9—30 

9  —  0.^3 

Half-Beam  Length  Ratio 

0.0S  —  0.03 

0.03  —  0.03 

Draft -Length  Ratio 

0.07  —  0.04 

0.0?  —  0.04 
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Thus,  when  we  consider  physical  quantities,  non-dimensionalized 
with  respect  to  L  and  V,  it  is  easy  to  see  that  ki  as  a  factor 
to  a  term  in  our  analysis  affects  the  order  of  magnitude  of 
the  term  in  the  case  of  ordinary  merchant  ships  since  it  is 
0(l/e).  On  the  other  hand,  for  fast  ships  kt  may  be  considered 
as  of  0(l ) . 

As  for  the  Green's  function,  we  may  separately  estimate 
the  magnitudes  of  the  local  disturbance  and  the  regular  part 
from  tne  following  simple  ship.  Namely,  when  we  consider  a 
wedge  ship  whose  waterline  slope  is 


HYDRONAUTICS,  Incorporated 


-28- 


which  is  0(1)  for  all  k*x> 


ki  2  -k;  d’ 

CPV  _J  *  2a  ^  min  e  ,  1  -  e  1 
x.reg!  o  / 


kx  x 

/v 


Y  (t)  at 
ov  7 


Since,  for  fast  ships,  kid  =  0(e),  we  o  tain  1  -  e 
this  case.  Therefore  we  can  writ0 


-ki  a 


=  0(e)  ir. 


cp 


x,reg 


=  C(e )  for  ’ow  and  medium  speeds 


^  «  0(ea  )  ior  nigh  speeds 


Thr  local  disturbance  on  the  free  surface  is  known  to  be  cf  tne 
same  order  rear  tr.e  booy  as  the  regular  wave  although  it  decays 
very  rapidly  In  the  far  field.  It  is  evident  from  the  express io 
for  +  that  the  derivative  with  respect  to  any  coordinate  non- 
oimenslor.alised  by  L  increases  the  factor  k>  while  the  integra¬ 
tion  reduces  this  factor. 

From  tre  first  twc  terms  of  Greer,  s  function  [ *3] , 


a  i- 

o  rx 


ra 


»  it  it:  not  difficult  to  show  that 


r»  ,~4  -i 

^  **  *  /  /  *«• 4X1 421 

o  i-C  c 


and  tp  are  Cue*  iog  f )  •  car  me  snip.  1*  is  also  easy  *u  see 
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that  the  continuous  variation  of  the  waterline  slope  does  not 
affect  the  order  of  magnitude  of  the  physical  variables  since 
the  influence  on  the  physical  quantities  from  the  sudden  change 
of  the  slope  at  the  bow  and  stern  dominates  the  other  terms  for 
ordinary  merchant  ships  and  is  of  the  same  order  as  the  other 
terms  in  the  case  of  fact  ships. 

From  the  boundary  condition,  we  know  that  cp^  is  of  the  same 

order  as  cp  /ki  on  the  free  surface  and  is  equal  to  zero  on  the 
xx 

ship's  flat  bottom.  We  also  know,  from  the  boundary  condition 
on  the  ship  surface,  that  cp  *  0(e )  near  the  ship  sides. 

y 

Thus  we  may  tabulate  the  magnitudes  of  some  physical  quan¬ 
tities  near  the  ship  surface,  except  within  a  distance  ea  from 
the  bow  and  stern  as  shown  in  Table  2. 

TABLE  2 


Order  of  Magnitudes  of  Some  Physical  Quantities 


Ordinary 
Merchant  Snip 

Fast  Ships 

1A-. 

0(c) 

0(1) 

CP 

C(ca  log  c  ) 

0(ca  log  €  ) 

v*2 

0(c) 

o(ca  ) 

Wave  Height  *  ®  /** 

A 

0(ca) 

o(«a  )  j 

cp 

V 

V 

0(c) 

0(€  )  1 
i 

V  .<p 

AA 

C{1) 

C(€*  )  j 

$  ,3 

yx  cx 

0(1) 

(C3  )  j 

E 
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Therefore  the  higher  order  potential  [28]  and  wave  re¬ 
sistance  [39]  in  the  case  of  slender  ships  are  only  valid  for 
fast  ships.  For  slender  ships,  at  low  Froude  numbers,  namely 
for  ordinary  merchant  ships,  the  higher  order  potentials  are 
rather  the  same  as  for  the  case  of  thin  ships  shown  in  [24]. 

The  lowest  order  non-dimensional  wave  resistance  for  or¬ 
dinary  merchant  ships  Is  the  same  as  that  given  in  Equation  [34], 
but  the  magnitude  should  be  listed  as  0(e4),  as  in  the  case  of 
fast  ships,  because  the  double  integral  in  Equation  [34]  brings 
out  the  factor  1/ki3  which  is  0(e2).  The  higher  order  wave  re¬ 
sistance  from  ordinary  merchant  ships  can  be  obtained  from  [32] . 
This  will  include  not  only  the  effects  needed  for  the  fast  ships 
in  [39]  but  also  the  other  effects  due  to  the  nonlinear  terms 
of  the  pressure  equation  [31]  viz  -I/2(cplx2  +  cp2  z2  )  which  is 
neglected  in  [39].  as  well  as  additional  terms  in  cp2x  shown  in 
l 24] .  Essentially  the  expression  for  the  higher  order  wave  re¬ 
sistance  of  ordinary  merchant  ships  will  be  the  same  as  for  the 
case  of  thin  snips.  However  if  kL  is  sufficiently  large  such 
that  e~klH  =  0(e),  (e.g.,  kx  =  30,  H  =  0.07,  then  e"klH  =  0.122) 
then  the  bottom  effect  and  the  effect  due  to  the  change  of  the 
wetted  surface  of  the  ship  in  [24]  will  become  one  order  higher. 

When  we  use  the  Taylor  series  expansion  of  Green : s  func¬ 
tion  [25],  we  have  to  notice  the  parameter  ki  e fi  where  the  factor 
ki  comes  In  because  of  a  differentiation  of  Green's  function 
wirh  respect  to  y.  The  parameter  kit  varies  in  the  range  of 
0,  fp  -  2.4  for  ordinary  merchant  ships.  Thus,  in  this  case, 

— e  convergence  of  [23]  is  unreliable  in  a  certain  part  of  fL  , 
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and  the  accuracy  of  the  first  order  theory  (Michell's  theory)  ap¬ 
pears  to  be  grossly  affected  and  the  convergence  of  the  higher 
order  theory  may  be  slow.  Nevertheless,  the  situation  is  not 
that  serious  due  to  the  fact  that  kiefi  in  [25]  is  a  function 
of  x  and  z  which  is  zero  at  x  =  0  and  x  =  1,  and  to  the  prin¬ 
ciple  that  the  potential  is  dominantly  influenced  by  the  entrance 
angle  for  low  Froude  numbers. 

SINGULARITIES  AT  THE  BOW  AND  STERN 

When  the  waterline  slopes  at  the  bow  and  the  stern  change 
suddenly,  cp^  has  a  logarithmic  singularity  at  the  origin  while 
cplx  and  cpiy  are  finite  there.  This  means  that  the  slope  of  the 
free  surface  and  F^  at  the  origin  is  infinity.  In  addition 

cplx  ana  the  higher  order  derivatives  are  singular  along  the  bcw 
and  the  stern.  Thus  the  present  analysis  cannot  hold  in  the  im¬ 
mediate  neighborhood  of  the  bow  and  stern.  However  the  princi¬ 
pal  effects  of  this  singularity  under  the  free  surface  is  the 
same  as  in  the  flow  about  a  thin  wing  In  an  infinite  medium  with¬ 
out  a  free  surface. 

This  fact  also  indicates  that  the  singular  behavior  of  the 
line  integral  [30a]  due  to  the  point  singularity  at  the  origin 
discussed  earlier  cannot  be  analyzed  by  means  of  this  linear 
theory.  In  fact  if  we  consider  the  line  integral  in  [23*1  or 
[30a],  in  the  physical  space  then  we  have 
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In  this  expression  any  effect  from  a  point  source  or  a  point 
doublet  on  the  free  surface  is  no  longer  included. 

CONCLUDING  REMARKS 

In  the  first  order  ship  theory,  whether  we  start  to  satisfy 
the  boundary  condition  on  the  free  surface  itself  or  on  the  mean 
free  surface  does  not  make  any  difference  in  the  result.  Like¬ 
wise,  whether  the  free  surface  is  considered  to  be  internally 
bounded  by  the  ship  waterline  at  the  free  surface  or  not  does 
not  affect  the  result  of  linear  theory.  However  when  we  develop 
a  higher  order  theory,  it  seems  necessary  to  consider  these  two 
items  very  carefully.  The  application  of  the  free  surface  con¬ 
dition  on  the  mean  free  surface  cannot  lead  to  the  proper  treat¬ 
ment  of  the  Influence  of  the  change  of  ship  wetted  surface 
caused  by  surface  waves.  This  difficulty  has  been  resolved  by 
distorting  the  coordinates.  Furthermore  the  requirement  that 
the  free  surface  be  internally  bounded  by  the  ship  waterline  at 
the  free  surface  makes  a  large  difference  in  the  result  obtained. 
The  aforementioned  influences  appear  in  the  second  and  higher 
order  terms,  and  are  more  Important  than  the  usually  considered 
nonlinear  effect  of  the  free  surface  boundary  condition  which 
comes  in  only  to  the  third  and  higher  order  terms  in  the  case  of 
fast  slender  ships,  and  to  the  second  and  the  higher  order  terms 
in  the  case  of  thin  ships  and  common  speed  slender  ships. 

The  higher  order  analysis  for  common  speed  slender  ships 
(ordinary  merchant  ships)  is  rather  the  same  as  that  for  thin 
ships,  and  the  equations  for  the  second  order  wave  height  and 
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for  the  second  order  wave  resistance  are  a  little  more  compile 
cated  than  those  for  fast  slender  ships. 

Essentially,  the  higher  order  wave  theory  considered  here, 
is  based  on  the  assumption  of  the  regular  or  asymptotic  con¬ 
vergence  of  the  series  expansion  of  physical  quantities.  In 
this  respect  any  higher  order  terms  can  be  obtained  from  the 
present  analysis.  However,  in  that  case,  unless  we  consider 
the  viscosity  and  the  surface  tension,  the  theory  higher  than 
the  second  order  would  be  meaningless. 

The  numerical  computation  of  the  higher  order  wave  form 
or  wave  resistance  due  to  a  slender  ship  is  not  too  complicated 
in  this  age  of  high  speed  computers.  It  would  be  extremely  in¬ 
teresting  to  see  the  numerical  results  even  for  the  case  of 
individual  high  order  influences  such  as  the  line  integral  ef¬ 
fect,  the  bottom  effect,  and  others. 
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